Quantum levitation of nanoparticles seen with ultracold neutrons 
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Analyzing new experiments with ultracold neutrons (UCNs) we show that physical adsorption 
of nanoparticles/nano-droplets, levitating in high-excited states in a deep and broad potential well 
formed by van der Waals/Casimir- Polder (vdW/CP) forces results in new effects on a cross-road of 
fundamental interactions, neutron, surface and nanoparticle physics. Accounting for the interaction 
of UCNs with nanoparticles explains a recently discovered intriguing small heating of UCNs in 
traps. It might be relevant to the striking conflict of the neutron lifetime experiments with smallest 
reported uncertainties by adding false effects there. 

PACS numbers: 



Surface diffusion of atoms, molecules and clusters, 
physically adsorbed in the potential well associated 
with van der Waals and Casimir-Polder interactions 
(vdW/CP) plays a key role in various phenomena in 
physics, chemistry, biology, and in applications [H-Q. 
The interaction is affected at small distance by "close- 
to-contact" effects depending on roughness and surface 
state [H|, Q . Because of experimental limitations the pre- 
cision studies are usually restricted to probes having ei- 
ther molecular or macroscopic sizes. Here we address 
the less extensively explored field of physical adsorption 
of nanoparticles and discover qualitatively new univer- 
sal phenomena related to their sizes and masses. Rigor- 
ous theoretical formalism describing states of physically 
adsorbed nanoparticles and their interaction with ultra- 
cold neutrons (UCN) 0^1 is given in a longer article 
(Ioj |; however general features can be understood from 
arguments stated below. The depth of the potential well 
affecting the motion of a sufficiently small physisorbed 
particle near a surface is proportional to the number of 
atoms in it, while the thermal energy attributed to the 
particle is always 3/2 k^T, where &b is the Boltzmann 
constant and T the temperature. Thus a nanoparticle 
is found near the surface in a deep and broad poten- 
tial well in normal direction; nanoparticles in states with 
low quantum number n are strongly bound to the sur- 
face, while those in high-n states move along it, with 
the effects of roughness and inhomogeneities mixing the 
two velocity components. Nanoparticles in high-n states 
might form a two-dimensional "cloud". The number of 
states is so large that a nanoparticle exhibits typically a 
quasi-classical motion fitting the conditions searched for 
in for solving puzzles in UCN physics. 

Here, we report on the theoretical justification for ex- 
istence of levitating nanoparticles, new treatment of the 
relevant existing UCN experimental data, and the results 
of new dedicated experiments, which provide evidence 
for levitating nanoparticles over solid and liquid surfaces. 



We show that our theoretical model describes all relevant 
existing and new UCN data; to our knowledge there are 
no alternative models, which could describe them. After 
a few reminders on the vdW/CP potential, we describe 
scattering of UCN on nanoparticles, discuss experimental 
results and consider consequences of our findings. 

The shape of the potential is calculated using a gen- 
eral expression involving only the scattering properties of 
the two objects [12j]. For example, the case of diamond 
nanospheres above a copper plane has been studied in 
great details in [l3| . Electromagnetic waves scattering on 
a plate is characterized by Fresnel reflection amplitudes, 
the form of which is fixed by the dielectric response func- 
tion s of copper. For frequencies uj lower than the plasma 
frequency o;p, e is large and copper behaves as a good 
reflector. For larger frequencies the reflection properties 
are poorer. When considering a nanoparticle of vanish- 
ing radius, the electric dipole approximation is sufficient. 
Calculations based on a simple model for the dielectric 
function of diamond containing only one resonance fre- 
quency uji reproduce the well-known vdW/CP energy. 
For distances small compared to the plasma wavelength 
(136nm for copper), or the resonance wavelength (106nm 
for diamond), the interaction reduces to the commonly 
used vdW formula with a power law R 3 /L 3 in the vicin- 
ity of the surface, where R is the nanoparticle radius, and 
L the distance of closest approach to the plane. 

For nanoparticles of arbitrary size the quantum disper- 
sion energy must be calculated by integrating the phase- 
shifts corresponding to all modes of the electromagnetic 
vacuum. In particular, many Mie scattering amplitudes 
(corresponding to many multipole components beyond 
the electric dipole), contribute significantly to the effect 
when the particle is close to the surface [l2j. The exact 
solution for the interaction energy predicts a smoother 
power law R/L close to contact with the surface, thus 
leading to a regular solution for the Schrodinger equa- 
tion, in contrast to the commonly used vdW formula. 
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We have analyzed the Schrodinger equation for the 
quantum dispersion potential and will sketch our find- 
ings here; details of the calculations are given in [Tq| . 
We found a large number of bound states with the 
eigenenergies E n deeper than 3/2/cb^, localized at dis- 
tances of up to a few nanometers to the surface, pro- 
vided R > lnm. The nanoparticle thus might fly over 
the surface while levitating in high-n quantum states. 
More precise statements depend on the population and 
thermalization dynamics and are sensitive to the low-n 
quantum states spectrum. In contrast to high-n states, 
where the nanoparticles stay essentially far from the sur- 
face, the precise properties of low-n states depend on 
close-to-contact details and in particular on the effect of 
roughness. Consequences of this effect can be charac- 
terized by using the methods developed for studying the 
Casimir force between two metallic plates. The probabil- 
ity of approach to close distances is reduced by contact 
repulsion from the highest peaks of the roughness profile 
[5[, while the dispersion force is affected by the distri- 
bution of approach distances due to roughness [6[. Here 
we chose a simple parametric approach: on top of the 
attraction calculated in [6[, we introduce contact repul- 
sion so that the effective potential has a minimum at 
a distance equal to l-3nm for the solid samples used. 
Because of the large number of high-n quantum states 
for the bound nanoparticles, a quasi-classical description 
of their motion in the vdW/CP potential is sufficient. 
Thus, when UCN bounce on the surface covered by lev- 
itating nanoparticles, we essentially deal with events of 
Doppler shift in UCN energy in the laboratory frame due 
to their elastic (in the centre of mass reference) collisions 
with nanoparticles in high-n quantum states. Due to ex- 
tremely low temperature of UCN (< lmK), small energy 
(< 10 _7 eV), and low velocity (a few m/s) UCNs possess 
the unique property of total elastic reflection from mo- 
tionless surfaces at any incidence angle. The probability 
of UCN loss per one bounce might be quite low: theoret- 
ically predicted probabilities of losses can reach ~ 10 -9 ; 
the best experimentally achieved values are ~ 10 -6 . The 
dominant loss mechanisms are nuclear absorbtion and 
up-scattering on phonons in the surface material. Even at 
the ambient temperature the probability of up-scattering 
is much smaller than the probability of coherent elastic 
reflection. 

That is why UCNs can be stored in closed traps for 
extended periods thus providing an extremely sensitive 
probe for rare processes or weak interactions. When 
studying UCNs storage in traps, unusual inelastic scat- 
tering of UCNs on trap surfaces was discovered 
One calls this phenomenon "small heating of UCNs" and 
such up-scattered neutrons Vaporizing UCNs (VUCN) - 
in analogy to vaporization of molecules. First measure- 
ments were followed by studies of several research groups, 
which essentially confirmed the initial observation, but 
suggested controversial estimations of the VUCNs pro- 



duction rates. 

Let us first give strong experimental motivation for 
the existence of levitating nanoparticles above surfaces. 
In the following we will show VUCN spectra that we cal- 
culated with a method sketched below using the initial 
UCN spectrum and the spectrometer spectral efficiency 
that we had measured experimentally. It has been found 
in [18| that the count rate of VUCN produced on stain- 
less steel surfaces as a function of the sample pre-heating 
temperature shows a sharp increase, the so-called "tem- 
perature resonance" shown in the left part of Fig.l. On 
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FIG. 1: (left) The probability of small heating of UCNs on 
surface of stainless steel samples is shown as a function of 
the temperature of sample outgasing (pre-heating); the mea- 
surement is performed at 300K. Black and white points show 
results of two independent experiments using analogous sam- 
ples, (right) The surface density of nanostructures with a 
radius of 6-7nm (stars) corresponding to the resonance en- 
hancement in the efficiency of VUCN detection, observed in 
AFM. More details in the body text. 

the right part of Fig.l we show the surface density of 
nanostructures as a function of temperature measured 
with an Atomic Force Microscope (AFM). These mea- 
surements indicate intense formation of nanostructures 
on the surface which takes place precisely at the temper- 
ature of sharp increase in the VUCN count rate. Nanos- 
tructures' size increases monotonously as a function of 
the heating temperature: a few particles at the ambient 
temperature, intense growth of the number and size of 
nanostructures below 350C; further increase in size but 
decrease in number of nanoparticles above 350C (some 
nanoparticles coagulate). While no alternative explana- 
tion has yet been proposed, we regard this as indication 
that it is due to levitating nanoparticles. 

Let us now briefly describe the basic equations of our 
model. The average differential cross-section of the in- 
teraction of an UCN with a nanoparticle is [Tq| : 

da _ 1 / M \ 3/2 
dE ~ 2mk n \2nk B T J * 

/ mv 2 \ r A;max r 27T 

dVeXP \2k^f)j k dko J MfB(k ,V,E^)\ 2 

with m, v n , k n = mv n the neutron mass, initial velocity 
and momentum and M, V the nanoparticle mass and 
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velocity. The kinematically allowed area of integration 
[^min^max] is determined by the conditions 



\mV — k n \ < ko < mV + k n 



\y / k 2 l + 2mE -mV\ < k < ^k 2 n + 2m E + mV (2) 

E is the energy transfer and ko = m\~& n — the incident 
momentum modulus in the center-of-mass reference sys- 
tem. The Born amplitude /b of the neutron nanoparticle 
scattering is /b(q) = —2mRUo/q 2 (hsm(qR/h)/(qR) — 
cos(qR/h)), where q is the transferred momentum, and 
Uo is the nanoparticle-neutron optical potential. To es- 
timate relevant nanoparticle radii, we use the above for- 
malism and follow the general approach from 1) R 
should not be too small, otherwise the coherent interac- 
tion cross-section is too low (as da/dE ~ i? 6 ), and 2) R 
should not be too large, otherwise nanoparticles are too 
heavy; thus their thermal velocity is too small, q is too 
small, and VUCNs could not be found in the window of 
efficient VUCN detection. 

We first apply our method to the small heating of 
UCNs on solid surfaces, where we have estimated in sev- 
eral manners the effective radius R of nanoparticles con- 
tributing to the observed inelastic process. The first sam- 
ple (la, LL7j we consider is powder of diamond nanopar- 
ticles with radii of 1-10 nm, on a copper surface, or 
above "sand" of diamond nanoparticles. We estimate 
the effective radius by the following three methods. 1) 
Using measured shapes of the differential VUCN spec- 

— spectrum * . 

tra leading to ^diamond = 9.4 ± 0.4nm. 2) Using mea- 
sured temperature dependencies of the count rates of 
detected VUCNs which gives ^diamond = 9.5 ± 0.6nm; 
3) Using the measured size distribution of nanoparticles: 

^diamond = 8.5 ± 0.5nm. Agreement between these inde- 
pendent estimations is convincing keeping in mind that 
they are not very sensitive to such parameters as the 
nanoparticle size and shape distribution, clustering, sur- 
face roughness, impurities etc. In contrast, if one aban- 
dons the "free levitation" model and "tunes" one critical 
parameter alone, say R by 10-20%, theoretical predic- 
tions largely contradict the data. Strictly speaking, there 
is no reason then, why all these estimations of R are 
equal, or even comparable. Another measured sample is 
nanoparticles appearin g on stainless steel surface due its 
thermal treatment [l3, LL8j. Again, the two independent 
estimations of the effective radii i^teei mm = 6.3 ± 0.3nm 
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and R stee \ = 6.6 ± 0.3nm agree. The effective masses 
of diamond and stainless steel nanoparticles estimated 
above are equal, thus providing an additional test of va- 
lidity of our model. 

Small heating of UCN on liquid Fomblin oils has been 
observed as well but no detailed spectral mea- 

surements have been performed. In view of the above 
evidence we assume that this phenomenon might be due 
to collisions of UCNs with levitating Fomblin oil nano- 
droplets. This mechanism would complement the known 



effects due to surface capillary waves [l9| and surface 
thermal fluctuations [lo[. We mention but not analyze 
here other hypotheses on origin of VUCN [2ll-[23| ; as they 
do not provide quantitative predictions and/or do conflict 
with experiments. While solid nanoparticles are charac- 
terized by long formation and evolution times, the num- 
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ber and size of levitating liquid nano-droplets 
could rapidly change with temperature, or other condi- 
tions like pressure, as they are governed by equilibrium of 
permanent formation of nano-droplets from vapors, also 
from explosion of bubbles on surface, and of their evap- 
oration. 

To compare these models, we measured VUCN count 
rates as a function of temperature, shown in Fig. 2. Al- 
though our model does not describe explicitly the num- 
ber of nano-droplets as a function of temperature, it is 
natural to assume that this number decreases with falling 
temperature as the Fomblin oil vapor pressure above sur- 
face, in analogy to the present data. Alternative inter- 
pretations of neutron small-heating, such as the hypoth- 
esis of thermal fluctuations of the surface (20j |. and the 
hypothesis of capillary surface waves [lflj, seem to con- 
tradict these experimental results, as they propose either 
too low absolute probabilities [20|, or /and inverse tem- 
perature dependence of the probability [l9j. In contrast 
our model gives the correct qualitative behavior. 
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FIG. 2: VUCN count rate in a.u. on a Fomblin surface as 
a function of temperature. The circle indicates the data at 
ambient temperature (the total probability of small heating 
is 10 -5 per wall collision); open squares show data measured 
when cooling the sample down, and open circles correspond 
to data measured when warming the sample up. 



In order to assess this question more precisely we have 
calculated the VUCN spectrum within our model and 
compared it to results of a new dedicated precision study 
of small heating of UCNs on solid and liquid surfaces 
in our precisely calibrated spectrometer BGS shown in 
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Fig. 3. The comparison is shown in Fig. 4 for several dif- 




FIG. 3: Layout of the Big Gravitational Spectrometer (BGS): 
(1) sample, (2) gravitational barrier, (3) entrance valve, (4) 
UCN monitor detector, (5) UCN absorber, (6) VUCN de- 
tector, and (7) exit valve. The principle of the procedure 
to measure small energy transfers is sketched in the insert on 
the left side. From bottom to top (in scale): initial differential 
UCN spectrum (the mean energy is 31neV , the half-width of 
the spectral mono- line is 3neV), a dead-zone of 3neV insen- 
sitive to VUCNs, the differential and integral (dashed line) 
spectra of VUCN. The differential efficiency of VUCN detec- 
tion is calculated using precisely measured values of UCN and 
VUCN storage times as a function of their energy; decrease in 
the detection efficiency, caused by partial losses of neutrons 
in samples, is taken into account. 

ferent physical systems, i.e. ultra-diamond (R ~2.5nm) 
and sapphire (R ~10nm) nanoparticles with broad size 
distribution, thin (1 jam) and thick (a few mm) layers 
of Fomblin oil, and naturally growing nanoparticles on 
a copper surface. Clearly all integral spectra are equiv- 
alent within statistical accuracy showing that there is a 
universal behavior for these very different physical sys- 
tems. This is a natural consequence of our model. In 
fact, the sensitivity of the spectrometer is sharply shaped 
to some nanoparticle mass, which, in turn, depends on 
the initial UCN energy range and the window of the spec- 
trometer sensitivity. However, within the peak resonance 
sensitivity, there is a slight dependence on the mass dis- 
tribution, which we show in Fig. 4. If the size distribution 
is known, the corresponding uncertainties are small. Ro- 
tations, shape and non-uniformity of nanoparticles, in- 
terference of scattering on neighbor nanoparticles and 
consequent scattering on one nanoparticle have been dis- 
regarded here. Concerning the results for Fomblin oil, 
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FIG. 4: VUCN count rate as a function of neutron energy, 
expressed in UCN raising height in the Earth's gravitational 
field, in cm. Circles show results measured with various solid 
nanoparticles: diamond, sapphire, copper; all these results 
agree within statistical accuracy. Rhombi indicate data mea- 
sured with a Fomblin oil sample; they agree with the data on 
solid nanoparticles. Four lines correspond to our model cal- 
culations of VUCN spectra on the Fomblin surface, for four 
hypotheses on the size distribution of nano-droplets: (i?/i?o) _ * 
where / = 1, 2, 3, 4 from top to bottom. The lines and the data 
are normalized to equal count rate at infinite height. 



they give evidence that the dominant mechanism of small 
heating of UCN on Fomblin oil surface in our experiment 
is UCN scattering on levitating nano-droplets. 

The present study was motivated in part by the unsat- 
isfactory status of the neutron lifetime r n experiments 
which show large and so far unexplained discrepancy be- 
tween results with smallest reported uncertainties 
l26j . This contradiction has been discussed with far- 
reaching consequences of an eventual shift of the mean 
world value for fundamental particle physics and cosmol- 
ogy (2?]-[29[, but only a single study has tried to v erify 
independently the validity of the measured results [30|. 
All these experiments use UCN traps with Fomblin-oil 
walls, and assume conservation of UCN energy. How- 
ever using our data and model, it is easy to show that 
a major fraction of UCNs do change their energy, so 
that their loss rate is different from the assumed val- 
ues and major false effects might arise. In particular, as 
the rate of forming Fomblin-oil nano-droplets depends on 
experimental conditions, the basic idea of these experi- 
ments is compromised: the geometrical and energy meth- 
ods of UCN loss extrapolation cannot be applied with- 
out reservations. Nano-droplet production rates depend 
on parameters never properly controlled like pre-history 
of the Fomblin oil treatment or vacuum. In particular, 
the UCN loss probability would be different in large and 
small traps; it would also evolve in time. The task of 
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estimating reliable corrections to r n values goes beyond 
the scope of this article. An attractive application might 
consist in studying vdW/CP interaction between levitat- 
ing nanoparticles and surfaces. Corresponding correc- 
tions to the integral VUCN spectra account for ~ 10% 
even in the present study, not optimized for these pur- 
poses. Surface potentials define the distribution of dis- 
tances of levitating nanoparticles to surface. As UCN 
could be up-scattered twice on one nanoparticle, before 
reflection from the surface and afterwards, the spectra 
would be modified depending on the distance to surface. 
Besides, complementary to using trapped atomic clusters 
for decorating surface defects, boundaries, step edges, 
grain boundaries, elastic strain fields, UCNs are sensitive 
to nanoparticles in motion above defect-free zones, thus 
giving us access to nanoparticle mobility. In levitation, 
nanoparticle mobility is very high; on the other hand 
chemical interactions largely reduce it thus giving us ac- 
cess to chemical properties of nanoparticles and surfaces. 
In an analogous way, one could also measure electrostatic 
effects. 

We thank our colleagues for useful discussions, in par- 
ticular A.L. Barabanov, S.T. Belyaev, H.G. Boerner, 
L.N. Bondarenko, D. Delande, D. Dubbers, R. Golub, 
V.K. Ignatovich, A. Leadbetter, S.K. Lamoreaux, V.I. 
Morozov, G.V. Nekhaev, S. Paul, J.M. Pendlebury, G. 
Pignol, K.V. Protasov. The authors benefited from 
exchange of ideas within the ESF Research Network 
CASIMIR, and the GRANIT collaboration. 



[1] K. Oura, V.G. Lifshits, A. A. Saranin, Surface science: 
An introduction. Bergin-Verlag, Berlin-Heidelberg : s.n., 
2003. 

[2] K.W. Kolasinski, Surface science: foundations of catal- 
ysis and nanoscience. Chichester, West Sussex : John 
Wiley & Sons, 2008. 

[3] G.Antczak, G. Ehrlich, Surf. Sci. Rep. 62 (2003) 39. 

[4] E.Shustorovich, Metal-surface reaction energetics: 



theory and applications to heterogeneous analysis, 
chemisorption, and surface diffusion, s.l. : VCN 
Publishers, Inc., 1991. 

[5] PC van Zwol et al, Phys. Rev. B. 80 (2009) 235401. 

[6] W. Broer et al, Europ. Phys. Lett. 95 (2011) 30001. 

[7] V.I. Luschikov, Yu.N. Pokoyilovsky, A.V. Strelkov, and 
F.L. Shapiro, JETP Lett. 9 (1969) 23. 

[8] V.K. Ignatovich, The Physics of Ultracold Neutrons. Ox- 
ford : Clarendon Press, 1990. 

[9] R. Golub, D.J. Richardson, and S.K. Lamoreux, Ultra- 
cold Neutrons. Bristol : Higler, 1991. 
[10] V.V. Nesvizhevsky, A.Yu. Voronin, A. Lambrecht, S. 
Reynaud, submitted to New J. Phys; arXivxond- 
matter.other/1205.1370 
[11] V.V. Nesvizhevsky, Phys. At. Nucl. 65 (2002) 400. 
[12] A. Canaguier-Durand et al, Phys. Rev. Lett. 102 (2009) 
230404. 

[13] A. Canaguier-Durand et al, Phys. Rev. A. 83 (2011) 
032508. 

[14] V.V. Nesvizhevsky et al, Europ. J. Appl. Phys. 6 (1999) 
151. 

[15] AV. Strelkov et al, Nucl. Instr. Meth. A. 440 (2000) 695. 
[16] L. Bondarenko et al, JETP Lett. 68 (1998) 691. 
[17] EV. Lychagin et al, Phys. At. Nucl. 65 (2002) 1995. 
[18] D.G. Kartashov et al, Int. J. Nanoscience. 6 (2007) 501. 
[19] S.K. Lamoreaux, R. Golub, Phys. Rev. C 66 (2002) 
044309. 

[20] A.L. Barabanov, S.T. Belyaev, Europ. Phys. J. B. 15 
(2000) 59. 

[21] A. Stepaniants et al, J. Low Temp. Phys. 113 (1999) 1159. 

[22] S.S. Malik et al, Phys. Lett. A 260 (1999) 328. 

[23] Y.N. Pokotilovki, Europ. Phys. J. B 8 (1999) 1; Phys. 

Lett. A 255 (1999) 173; Y.N. Pokotilovki, JETP Lett. 69 

(1999) 91. 

[24] S. Arzumanov et al, Phys. Lett. B. 483 (2000) 15. 
[25] A.P Serebrov et al, Phys. Lett. B. 605 (2005) 72. 
[26] A. Pichlmaier et al, Phys. Lett. B. 693 (2010) 221. 
[27] H. Abele et al, Phys. Rev. Lett. 88 (2002) 211801. 
[28] D. Dubbers, M.G. Schmidt, Rev. Mod. Phys. 83 (2011) 
1111. 

[29] G.J. Mathews et al, Phys. Rev. D. 71 (2005) 021302. 

[30] S. Paul, Nucl. Instrum. Meth. 611 (2009) 157. 

[31] K. Koga et al, J. Chem. Phys. 109 (1998) 4063. 

[32] T. Tadros et al, Adv. Col. Interface Sci. 108 (2004) 108. 



